Universal macroscopic background formation in surface super-roughening 
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We study a class of super-rough growth models whose structure factor satisfies the Family- Vicsek 
scaling. We demonstrate that a macroscopic background spontaneously develops in the local surface 
profile, which dominates the scaling of the local surface width and the height-difference. The shape 
of the macroscopic background takes a form of a finite-order polynomial whose order is decided from 
the value of the global roughness exponent. Once the macroscopic background is subtracted, the 
width of the resulting local surface profile satisfies the Family- Vicsek scaling. We show that this 
feature is universal to all super-rough growth models, and we also discuss the difference between the 
macroscopic background formation and the pattern formation in other models. 
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Recently, there has been considerable progress in the understanding of the dynamics of growing surfaces Much 
of the progress is motivated by the crucial observation that surface roughening exhibits scaling behaviors. For example, 
the growth of the global surface width w{L,t), when it starts from the flat surface, shows the following behavior. 
At an initial stage, w{L, t) grows as a power law of time w{L, t) ~ , where fi is called the growth exponent. At a 
later stage, however, w{L,t) saturates to a certain power of the system size L, w{L,t) ~ L", where a is called the 
roughness exponent. The crossover to the saturated surface is governed by the lateral correlation length ^{t), which 
scales as t^/^ for the initial stage t <^ and saturates to L for the later stage t ^ L^. Here the third exponent z is 
called the dynamic exponent, and this scaling behavior is called the Family- Vicsek (FV) scaling ansatz 

The presence of the scaling is usually related to the self-affine structure of a surface, implying the same properties 
at all length scales. In some growth models (for example, Refs. [QH]), however, surface properties exhibit different 
behaviors when they are probed at different length scales. For example, one may probe the surface width w{l,t) 
within a local window of size I L), where 

w'il,t) ^ {{{h{x,t) ~ {hix,t))x, f)x,} , (1) 

and ( )xi is the spatial average over the local window Xi of size / and { } the sample average. And one finds 

r for t < P 

w^i,t)ocl , fmV^""""' ■ (2) 



Note that while it follows the FV scaling at the initial stage t <^ F, the scaling behavior at the later stage is anomalous 
: the window size dependence is described not by the global roughness exponent a, but rather by a new exponent 
a'{< a) which is called local roughness exponent. Note also that at the later stage, the surface width implied by 
Eq. is anomalously larger than the FV scaling result P", since a > a' and (,{t) ^ I. 

In the past few years, there have been many works on the anomalous roughening. Numerical simulations demon- 
strated the occurrence of the anomalous roughening in some growth models (for example, see Refs. QJ^). New scaling 
ansatz for the anomalous roughening was proposed . Very recently, Lopez et al. identified two separate mech- 
anisms of the anomalous roughening, super-roughening and intrinsic roughening. In case of the super-roughening 
(a > 1), the structure factor S{k), or the power spectrum, follows the Family- Vicsek scaling and the anomalous 
roughening occurs due to the divergence of S{k) near fc = 0. In case of the intrinsic anomalous roughening, on the 
other hand, S{k) itself scales anomalously. 

Among the two mechanisms, we focus in this paper on the anomalous roughening due to the super-roughening. In 
particular, we study the morphology of the local surface profile generated in the super-rough growth models. Usually, 
the surface profile can be probed from the scaling of the local surface width or the height-difference as a function 
of a probing length scale. In super-rough surfaces, however, the scaling is much less informative, since the relevant 
local roughness exponent a' is always 1 Q independent of many details of growth models, the only implication of the 
universal value 1 being the divergence of S{k). So to probe the surface profile in the super-rough growth models, we 
take a different theoretical approach based on the least-square-fitting method. As a main result of the paper, we find 
that (i) the local surface profile forms a macroscopic background that takes the form of a finite order polynomial and 
that (ii) the magnitude of the short wavelength fluctuations superimposed on the macroscopic background follows the 
FV scaling. We also show that this property is universal to all super-rough growth models. 
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One simple way to study the profile is to make a guess on the functional form of the profile, fit the function to the 
profile through the least-square-fitting method, and examine the magnitude of the discrepancy. The simplest example 
of this approach is the local surface width since the idea of Eq. (|l|) is nothing but the least-square-fitting of the local 
surface profile with a constant function {h(x,t)) x, ■ Adopting this view point, it is then interesting to explore other 
possibilities of fitting functions, which hopefully produces a smaller discrepancy. Since the origin of the anomalous 
roughening in the super-rough growth models is the divergence in the long wavelength components that appear smooth 
in the local window, we first take a first-order polynomial in x as a fitting function : hi{x,t) = ai(t)x + ao{t). Here, 
the coefficients ai{t) are fixed by minimizing the discrepancy {[h{x,t) — hi{x,t)]'^)xi for a given surface height profile 
h{x, t). 

For illustration, we first use the 1+1 dimensional uniform diffusion models as examples : 

+ (3) 

where 77(0;, t) is the noise with correlation t)ri{x' , t')} =D5{x — x')S{t — t'), and m = 1, 2, 3, . . . . For an arbitrary 
m, it can be explicitly verified that the structure factor follows the FV scaling. Also the precise values of the scaling 
exponents are known : a = {2m — l)/2,/3 = {2m — l)/4m, and z = 2m Notice that for m > 2, the uniform 
diffusion models become super-rough (a > 1). Then as demonstrated in Ref. the models exhibit the anomalous 
scaling with a' = 1 for m > 2. 

The accuracy of the fitting with the first-order polynomial can be estimated from the ensemble averaged minimal 
discrepancy, 

wl{l,t) ^ {{{h{x,t) - hi{x,t)f)x,} , (4) 

where we call wi{l,t) first-order generalized surface width. For the linear models Eq. (||), wi{l,t) can be evaluated 
analytically in a straightforward way. At the initial stage (,{t) ^ /, w'f{l,t) scales the same way as w'^{l,t), implying 
that there is no substantial improvement of the fitting from the new choice of the fitting function. This result is 
understandable since the divergence of S{k) does not qualitatively affect the surface profile at the initial stage. At 
the later stage ^{t) ^ /, however, the results are 

r for a < 2 (or m = 1,2) 

wl{l,t) oc I , fmV^"^^^ , ^ (5) 



\ for a > 2 (or m > 3). 

Here the reduction in the power of the large factor S,{t)/l should be noticed : it has been reduced from 2(a — 1) in 
w'^{l,t) to for m — 2 and to 2{a — 2) for m > 3. For the definiteness, let us take l/£,{t) as a small parameter of 
the analysis. Then the power reduction implies the substantial improvement of the fitting, which we interpret as an 
evidence of a spontaneous formation of a macroscopic background in the surface profile. Here the word "macroscopic" 
denotes that the characteristic length scale of the background is comparable to or larger than the window size 
Figure |l| shows a typical surface profile of the m = 2 uniform diffusion model at the later stage (,{t) ^ I, and the 
formation of the linear background is clear ||l0[| . 

Motivated by the success of the first-order polynomial, we explore the idea further and examine the other types 
of fitting functions. Specifically we choose a higher-order polynomial, which is a natural extension of the first-order 
polynomial. One merit of this choice is that polynomials form a complete set of basis for functional space. So by 
examining the fitting with high-order polynomials, one can examine the effect of all types of fitting functions. To 
estimate the accuracy of the fitting with an A^th order polynomial, we introduce an A^th order generalized surface 
width WN{l,t), where 

w%{l,t) = {{[h{x,t) - hN{x,t)]^)x,} , (6) 

and the coefficients ai{t) of h]\r{x,t) = X^i^o '^ii£l^^ similarly fixed by the least-square-fitting method for a given 

eal 

w'^{l,t) at the initial stage, again finding no improvement, and 



realization h{x,t). For the linear model Eq. (|^), one can analytically verify that wj^{l,t) scales the same way as 



for a < A'' -I- 1 (or m < A' -I- 1) 
A^V ' ^ ^ / ili J for a > A^ + 1 (or m > A^ + 2) ' ^' 
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at the later stage. This scaUng result contains very interesting information. First, the accuracy improves substantially 
as the order N increases till N becomes larger than a — 1, after which it saturates and does not improve upon the 
further increase of N. Recalling the completeness of the polynomial basis, one then realizes that even if other forms of 
fitting function were chosen, such as sines and cosines, they would not improve the accuracy either. In this sense, one 
can say that the macroscopic background in the uniform diffusion models takes the form of an iVth order polynomial, 
where N is the largest integer smaller than a. Second, the minimal discrepancy of the fitting is always of the order 
of 0, and so it is characterized by the global roughness exponent. One interesting consequence of this result is 
that once the macroscopic background is subtracted, the width of the resulting local surface profile satisfies the FV 
scaling. 

To understand the underlying physics in the uniform diffusion models, let us represent the local surface profile in 
the following way : h{x,t) = Y^^=o ^ '^ii^)^^ + t) sin /ca; + i?(fc, t) cos fcx], where kl is an integer multiple of 

2tt. One interesting property of the uniform diffusion models is that while the relaxation force tends to suppress the 
amplitude of the sinusoidal part, it does not affect the polynomial part : F[Y^'^^f^^ ai{t)x'^] = for arbitrary ai{t) 
where F[h{x, t)] = (— This property opens a possibility for the anomalous growth of the coefficients 
ai{t), which makes polynomials special candidates for a fitting function. And the scaling result of w'j^{l,t) implies 
that the coefficients ai{t) with i < a (or i < m) indeed become anomalously large so that they lead to the formation 
of the macroscopic background. 

Let us compare the formation of the macroscopic background in the uniform diffusion models with pattern for- 
mation in other models where the pattern formation is due to symmetry breaking or the coherent superposition of 
modes of different wavelength In the uniform diffusion models, the symmetry breaking does not occur and the 
ensemble average of ai{t) vanishes. Also there is no coherent coupling between different modes since the models are 
linear. Here, however, the fluctuations of ai(t) in different realizations are so large that the fluctuations lead to the 
macroscopic background formation for any given realization, that is, the macroscopic background formation is due to 
large fluctuations in some particular "degrees of freedom" . In this sense, the macroscopic background formation is 
completely different from the pattern formation in its origin. 

For the 1+1 dimensional uniform diffusion models, it is now demonstrated that (i) a macroscopic background 
develops spontaneously in the local surface profile, and the background takes a form of an A^th order polynomial 
where N is the largest integer smaller than a, and that (ii) once the macroscopic background is subtracted, the 
resulting local surface profile satisfies the FV scaling. Below, we show that these properties are universal to all 
super-rough growth models. 

From now on, we restrict the discussion to the later stage since it is trivial to show that w%{l,t) satisfies the FV 
scaling at the initial stage for an arbitrary A^. To study the scaling behavior at the later stage, it is useful to first 
relate the height-height correlation function G{x,t) to the structure factor S{k,t) — {h{k,t)h{~k,t)}. By comparing 
the definitions, one can obtain the following relation between G{x,t) and S{k,t) ^ : 

G{x,t)=A —[\-cos{kx)]S{k,t), (8) 
ho 27r 

where the lower cutoff fco of the integration is decided by the system size L, fco — 2-k/L. Usually the lower cutoff 
can be replaced by and Eq. (||) leads to G{x,t) satisfying the FV scahng. For a > 1, however, the integration has 
the infrared divergence and the lower cutoff becomes important. In this case, by taking explicit care of the lower 
cutoff and from the FV scaling of S{k,t) [for the FV scaling of S{k,t), see, for example, Ref. Q], one can obtain for 

G{x,t) = |a;p"M(x/ii/") + \x\''^^'^''-^\t)v{x/^{t)), (9) 

where u{y 0) = uq [uq log(l/i/) if a is an integer], and v{y) is an even function of y, v{y) = X^i^o ^2i2/^*- In Eq. (|^), 
the first term is the usual Family- Vicsek contribution and the second is the anomalous contribution coming from the 
integration near the lower cutoff. The even form of the scaling function v{y) originates from the series expansion of 
the only a;-dependent factor 1 — cos(fcx) in Eq. (||). Notice that the anomalous contribution is dominant over the FV 
contribution for a > 1. 

To demonstrate the universality, we show that the result in Eq. , regarding the scaling of wf^ (l,t), comes naturally 
from the form of G{x, t) in Eq. (||). For simplicity, let us take A^ = 1. The generalization to a larger value of A^ is also 
straightforward and it will be sketched below. The key element of the demonstration is to derive a general relation 
between w\(l,t) and G{x,t). From the definition of w\{l,t), one can verify the following linear relation : 

!/■,.„. . f 12 



w^{l,t) = ^j dxidx2 G{xi - X2,t)\l + -1^x1X2] , (10) 
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where the center of the local window Xi is chosen as the origin of the coordinate system. Then by combining Eq. ( |lO| ) 
with Eq. (||), one finds 



i=0 



2(a-l-i) 



where 



^'^^^^ (2/i+l)(2M + 2)(2Ai + 4) ' ^^^^ 

and the dots represent other terms smaller than P". Then to find out the scaling behavior of w'f{l,t), one just needs 
to select the leading order term from Eq. (11). In power counting of l/^{t), /i(l)wo^^^(i)^*'"~^^ appears to be the 
leading order term. However this term is not the true leading order term, since the proportionality constant /i(l) 
vanishes identically. Then by choosing the next leading order term in power counting, one can show that Eq. (|^), 
which is a special case of Eq. (^ , is universal. 

The generalization of the analysis to a larger value of N can be carried out in a similar way. From the definition of 
w'j^{l,t), one can derive a linear relation between wff(l,t) and G{x,t). By combining this relation with Eq. (^, one 
finds 



"'jv(^ = 2^ /jv(* + l)«2^r" I ^ 1 + InHuoI'" + ■■■ , (13) 

where the dots represent other terms smaller than and /Ar(//) is zero if /z is an integer smaller than or equal to 
N |l3|| . Then by selecting the true (nonvanishing) leading order term from Eq. ([l3|), one can show the universality of 
Eq. (I 

So far, we have demonstrated the universality of the macroscopic background formation only for 1+1 dimensional 
systems. However the universality is not restricted to 1+1 dimensional systems. For illustration, here we show the 
universality of Eq. (||) for D+1 dimensional systems. The generalization of the two key relations Eqs. (|^,p^ can be 
achieved through trivial replacements such as fca; — > k • x and X1X2 Xi • X2, where a D-dimensional cube with 
volume is chosen as a D-dimensional local window. Then following the same procedure, one obtains precisely the 
same expression as Eq. ( [ll| ) except for a replacement of fi{fJ.) by /i,d(a*) = Dfi{fi), and one verifies the universality 
in higher spatial dimensions. The generalization of the analysis to a larger value of N is also straightforward. 

One important signature of the super-rough growth models is the universal value of the local roughness exponent 
a' ~\. For some growth models numerical simulations resulted in a' « 1, and so it would be interesting to 

examine whether the present analysis applies to those growth models. One of the most promising candidates is the 
1+1 dimensional model of driven interfaces in random media, for which the functional renormalization group study 
by Narayan and Fisher resulted in a' = 1 in the critical region. 

In this paper, we consider the anomalous roughening in super-rough (a > 1) growth models whose structure factors 
satisfy the Family- Vicsek scaling. We show that these growth models share the following universal features. First, 
the local surface profile is characterized by a formation of a macroscopic background, which takes a form of an A^th 
order polynomial, where N is the largest integer smaller than the global exponent a. Second, after subtracting the 
macroscopic background, the width of the resulting local surface profile satisfies the Family- Vicsek scaling. 
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FIG. 1. Typical surface profile of a 1+1 dimensional uniform diffusion model with m — 2 [Eq. (|3|)] at the later stage ^(t) 3> I- 
Notice the formation of a linear macroscopic background hi{x,t) whose large slope leads to the anomalous scaling of w^{l,t). 
Once the background is subtracted, the resulting local surface width wi{l,t) is characterized by the global roughness exponent 
a and satisfies the Family- Vicsek scaling. 
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